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Abstract. Trap models are intuitively appealing and often solvable models of glassy 
dynamics. In particular, they have been used to study aging and the resulting 
out-of-equilibrium fluctuation-dissipation relations between correlations and response 
functions. In this note I show briefly that one such relation, first given by Bouchaud 
and Dean, is valid for a general class of mean-field trap models: it relies only on the way 
a perturbation affects the transition rates, but is independent of the distribution of trap 
depths and the form of the unperturbed transition rates, and holds for all observables 
that are uncorrelated with the energy. The model with Glauber dynamics and an 
exponential distribution of trap depths, as considered by Barrat and Mezard, does not 
fall into this class if the perturbation is introduced in the natural way by shifting all trap 
energies. I show that a similar relation between response and correlation nevertheless 
holds for the out-of-equilibrium dynamics at low temperatures. The results points to 
intriguing parallels between trap models with energetic and entropic barriers. 



1. Introduction 

Trap models consist of a single particle, or equivalently an ensemble of non-interacting 
particles, hopping in a landscape of traps of energy E. Such models have been studied 
extensively and shown to account qualitatively for many interesting features of glassy 
dynamics, see e.g. PQI2E1I110EIIZ1IHIE1 CD]- In the simplest case the rates for 
transitions from one trap to another depend only on the energies of the two traps. One 
then has a mean-field trap model, where no information on any spatial organization is 
retained. This is the case that will concern us here; for work on spatial trap models see 

e.g. nu na nu - 

In this paper I focus on the behaviour of trap models after a quench into the 
glassy phase, and in particular on two-time correlation and response functions. For a 
generic observable m the correlation function is C(t, t w ) = (m(t)m(t w )) — (m(t)) (m(t w )), 
while the (linear) response function x(t,t w ) measures the change in (m(t)) due to a 
conjugate field h that is switched on at time t w < t. The time of preparation of the 
system by quenching is taken as the zero of the time axis, so that the waiting time 
t w can alternatively be thought of as the "age" of the system at the time when the 
field is applied. Over the last decade it has been recognized that out-of-equilibrium 
fluctuation-dissipation (FD) relations between such correlation and response functions 
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are extremely useful for characterizing glassy dynamics [To"! ITE! IT7| ITH] . One defines a 
fluctuation-dissipation theorem (FDT) "violation factor" X(t,t w ) by 

9 -x(Mw) = ^%^^C(M w ) (1) 



dt w ^ ' wy T <9t ¥ 
so that X = 1 corresponds to the usual equilibrium FDT. The value of X can be read 
off from the slope of a parametric "FD plot" of x versus C; see [HI EI] for a discussion of 
the effects of using either t or t w as curve parameter. The quantity on the l.h.s. of (JTJ), 
denoted R(t,t w ) below, is the response to a short field impulse at time t w . 

In glassy systems one typically finds that the decay of two-time correlation functions 
C(t,t w ) exhibits several regimes: an initial decay to a plateau, with further relaxation 
taking place only on "aging" timescales t — t w that grow with the age t w , for example 
t — t w = 0{t w ). In mean-field spin glasses [T51 HEl HZj one finds that, in the limit 
of large t w , X has well-defined and distinct values in these regimes, corresponding to 
an FD plot made up of two straight line segments: X = 1 in the short-time regime, 
corresponding to quasi-equilibrium, and X < 1 in the aging regime. In the latter, one 
can then define an effective temperature by T e g = T/X. This has been shown to have 
many of the properties of a thermodynamic temperature [151 CHI HZj, opening up the 
exciting prospect of an effective equilibrium description of out-of-equilibrium dynamics. 

The existence of effective temperatures in systems other than the now canonical 
(e.g. spherical p-) spin glass models has been the subject of much research in recent 
years, but a coherent picture has yet to emerge [IB]. In Bouchaud's trap model [I], 
intriguing results have recently been found 0: even though the correlation functions 
C(t,t w ) decay within a single aging "time sector" t — t w = (D(t w ), the FDT violation 
factor X is not constant as one might expect by analogy with mean-field spin glass 
models. Instead it varies continuously with (t — t w )/t w , resulting in a curved FD plot 
with an asymptotic slope (for (t — t w )/t w — > oo, i.e. C — > 0) of X^ = 0. 

In the Bouchaud trap model glassy dynamics arises from the presence of energy 
barriers. The aim of this paper is to analyse the FD relations in a different trap model, 
due to Barrat and Mezard 0, where glassiness instead results from the presence of 
entropic barriers. Because of the different physical mechanisms causing the out-of- 
equilibrium behaviour, it is then not a priori clear whether the FD relations of the two 
models should be related. I find that some important aspects of the FD relations are 
indeed the same, pointing to intriguing parallels between models with energetic and 
entropic barriers that deserve to be explored further. 

Very recently, Ritort ^0) has also considered FD relations in the Barrat and 
Mezard model. However, he assumed that the effect of the perturbing field on the 
transition rates defining the trap model dynamics has a simple multiplicative form. 
This is "not easy to justify a priori 1 JU] and gives only an approximation to the 
natural prescription where the effect of the perturbing field is to shift all energies 
according to E — > E — hm. I show in this paper that the response in this natural 
Barrat and Mezard model can be analysed directly, and I give exact results for the 
FD relations in the limit of low temperatures; these differ in important respects from 
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those obtained with the approximation of multiplicatively perturbed rates. As a by- 
product of the calculation, I also show briefly that the FD relation due to Bouchaud 
and Dean [3] , which was recovered by Ritort [10] for the Barrat and Mezard model with 
the approximation of multiplicatively perturbed rates, is in fact valid for arbitrary trap 
models with multiplicatively perturbed rates. 

In the following section I give the definitions of the Bouchaud and the Barrat and 
Mezard trap models; general expressions for correlation and response functions that 
apply to all trap models are then derived in Sec. El Sec. |U applies these to the case of 
multiplicatively perturbed rates. Sec. El contains the main result, the exact low-T FD 
relation for the Barrat and Mezard model. I conclude in Sec. IHlwith a discussion of the 
intriguing links between trap models with energetic and entropic barriers which arise as 
a consequence of this relation. 

2. Trap models 

A trap model is defined by a distribution p(E) of trap energies; the convention for the 
sign of E is here that lower E corresponds to deeper traps, which is the reverse of that 
in e.g. [3]. The primary dynamical quantity is then Po(E,t), the distribution of finding 
the particle in a trap of energy E at time t; the subscript indicates that for now we are 
considering the dynamics without any perturbing fields. The evolution of Pq is given by 
the master equation 

-P (E,t) = -T {E)P {E,t)+p{E) JdE'w (E <- E')P {E',t) (2) 

where wq(E <— E') is the rate for transitions between traps of energy E' and E. More 
precisely, if one considers a finite number of traps N, the transition rate from trap i 
to j is (l/N)wo(Ej <— Ei); the total rate for transitions to traps in the energy range 
E < Ej < E + dE is then Wq(E <— Ei) times the fraction of traps in this range, which 
is p(E)dE for large N. The quantity 

T (E) = J dE'p(E')w (E' <- E) (3) 

in (J2J is the total "exit rate" from a trap of energy E. 

Two specific instances of trap models have received considerable attention in recent 
years. Bouchaud pQ chose for his trap model p(E) = T g 1 exp(E/T g ) with — oo < E < 0. 
For any choice of transition rates that satisfies detailed balance, the model then has a 
glass transition at T = T g since the equilibrium distribution P eq (E) oc p(E) exp(—E/T) 
becomes unnormalizable there. For lower T, the system must show aging, i.e. a strong 
dependence of its properties on the waiting time t w . Bouchaud pQ| assumed transition 
rates 



w (E' <- E) = exp(pE) 



(4) 



Fluctuation- dissipation relations in trap models 



4 



that are independent of the energy of the arrival trap; here (3 — 1/T as usual. Barrat 
and Mezard |2] chose instead Glauber rates 

As emphasized by Ritort ^UJ > the out-of-equilibrium dynamics of these models is rather 
different: in the Bouchaud model with its activated dynamics, glassiness arises from the 
presence of energy barriers, and the system arrests completely for T — > 0. In the Barrat 
and Mezard case, on the other hand, the system can keep evolving by transitions to 
traps with ever lower energies, even at T = 0; the diminishing number of such traps 
effectively creates entropic barriers that slow the relaxation. 

3. Correlation and response 

We now want to consider the correlation and response properties of some, essentially 
arbitrary, observable m. In the most general terms the properties of this are described 
by the distributions p(m\E) of m across traps of given E. I will assume throughout that 
m is on average uncorrelated with E, so that its conditional mean 



= Jdmmp(m\E) (6) 
vanishes for all E; the variance 

A 2 (E)= fdmm 2 p(m\E) (7) 



however can be dependent on E. With m included, the master equation is 
^P(E,m,t) = -T(E,m)P(E,m,t) 

+ p{m\E)p(E) J dE 1 dm 1 w{E, m <- E', m')P(E', m', t) (8) 

where the rates w(E,m <— E',m') may now depend on a perturbing field h conjugate 
to m, and the total exit rates are 

T(E,m) = jdE' dm' p(m'\E')p(E')w(E',m' <- E,m) (9) 

An expression for the correlation function of m is easily found. In the absence of 
a field, w(E',m' <— E,m) = Wo(E' <— E) and T(E,m) = Tq(E) are independent of the 
value of our observable. Equation (jHJ) then shows that P(E,m,t) = p(m\E)P (E,t) as 
long as the same is true at time t = 0. (This is a natural assumption and holds e.g. 
when P(E, m, 0) is an equilibrium distribution at zero field and some initial temperature 
above T g , from which the system is quenched to T < T g at t = 0.) For our zero mean 
observables (JHJ) this implies in particular that (m(t)) = at all times. The two-time 
correlator of m is then 

C(t,t w ) = (m(t)m(t w )) (10) 

= JdEdmdE' dm'mm'P (E 7 m\E',m',t-t w )p(m'\E')P (E',t w ) (11) 
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Here P (E,m\E',m',t — t w ) is the propagator, i.e. the probability of being in a trap 
with energy E and observable m when starting from a trap with E' and w! a time t — t w 
earlier. This can be obtained as the solution to (JHJ) starting from the initial condition 
8{E — E')5(m — m'). Since the correlation function is calculated in the absence of a 
field, the only nontrivial m-dependence in (JHJ) arises from the factor p(m\E). Treating 
the second term on the r.h.s. of (JHJ) as an inhomogeneity one thus sees that 

P (E, m\E', rri, t - t w ) = e - r ° {E ' ){t - Uv) 5(E - E')8{m - m') + p(m\E) x ... (12) 

where the dots indicate factors not involving m. Inserting into (jllj) and using the 
zero-mean assumption (JHJ) then yields the simple representation 

C(t,t w )= JdE A 2 (E)e- r °^ t - t ^P (E,t w ) (13) 

for the correlation function. This makes sense: physically, every hop completely 
decorrelates the observable, so that C is an average of the probabilities exp[—T (E)(t — 
t w )} of remaining in the current trap, weighted by the probability of being in a trap of 
energy E at time t w and multiplied by the variance of m across traps of this energy. 

To find the impulse response _R(t, t w ), consider a field impulse of amplitude h and 
infinitesimal length At , applied at time t w . Denote 

Aw(E', m <- E, m) = w(E', m! *-E,m)- w (E' <- E) (14) 

the change in the transition rates caused by the field, and AT(E, m) similarly the change 
in the total exit rates; /i-dependences are not written explicitly here. Then from the 
master equation (JHJ), and using that P(E,m,t w ) = p(m\E)P Q (E,t w ), one has 

P(E, m,t w + At)= p(m\E)P (E, t w ) 

-AtT (E)p(m\E)P (E,t w ) 

+ Atp(m\E)p(E) JdE'dm'w (E <- E')p(m'\E')P Q (E' ,t v ) (15) 
- At AT(E, m)p{m\E)P (E, t w ) 

+ At p(m\E)p{E) J dE' dm' Aw{E, m «- E', m')p(m'\E')P (E f , t w ) 

where the effects of the field have been explicitly separated off in the last two lines. 
After time t w + At , when the field is switched off again, the same argument that lead 
to (JHJ) applies and so 

P(E, m, t) = e - Vo{E){t -^ ) P{E, m, t w + At ) + p{m\E) x ... (16) 

for t > t w + At with the dots again indicating factors independent of m; in the exponent 
I have approximated t — t w — At w t — t w since we are interested in the limit At — ► 0. 
To find (m(t)), from which the response function is obtained, one inserts ()15j) into (jlfij) . 
multiplies by m and integrates over m and E. All terms of the form p{m\E) x . . . give 
a vanishing contribution due to (JHJ). Only the last two lines of (JTKJ) thus survive, and 
the two-time response function can be written as 

hR(t,t w ) =±^( m (t)) (17) 
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= JdEdme-^^-^m -AT(E,m)p(m\E)P (E,t w ) 

+ p{m\E)p{E) JdE' dm' Aw(E,m <- E' ,m')p(m'\E')P (E' ,t v 

So far this applies for arbitrary field amplitude h, so that R(t, t w ) is in general a nonlinear 
response function, but we will specialize to the linear response limit h — > below. 

4. Multiplicatively perturbed rates 

To get concrete expressions for the response function one needs to define how the field 
h affects the transition rates. The natural prescription is that all energies are shifted 
according to the value of the observable, E — > E — hm and E f —> E' — hm'. Before 
going on to consider the more complicated case of the Barrat and Mezard model, I first 
briefly review the situation in the Bouchaud model, where a simple relation between 
correlation and response exists JI]. The derivation will show that this relation actually 
applies rather generally, being dependent only on the way the field affects the transition 
rates. 

For the Bouchaud model, shifting the energy E — > E — hm in (jlj) gives the 
transition rate in the presence of a field w(E',m' <— E,m) = exp(/3E — f3hm) = 
exp(—/3hm)wo(E' <— E). More generally, one can consider rates perturbed by the 
field according to [I] 

w{E\ m' <-E,m) = e ph[{1 - c)m '~ Cm] w (E' <- E) (19) 

which reduces to the natural choice^ for ( = 1 but also maintains detailed balance for 
other values of £. To linear order in h one then has 

Aw(E', m! *-E,m) = (3h[{l - Qm' - (m}w (E' 4- E) (20) 

and the corresponding change in the exit rates is 

AT(E,m) = (3h I dE' dm' [(1 - Qm' - (m]p(m'\E')p(E')w (E' <- E)(21) 

Using again the zero mean assumption (|6j) . the first term in square brackets vanishes, 
giving with (JHJ) 

AT(E,m) = -{3h(mT (E) (22) 

One can now substitute (|2T)|) - with the arguments (E, m) and (£", m') interchanged 
appropriately - and (}2*2|) into (fTHj) . Using fdm'm'p(m'\E') = and dividing by h yields 
for the linear response function 

-r {E)(t-t w ) 



R(t,t w ) = (3 JdEA 2 (E)e~ 



(T (E)P (E,t v 



+ (1 



C)p(E) JdE'w (E^E')P (E',t v 



(23) 



I The multiplicative perturbation of rates (|19|l arises from the natural energy shift prescription 
E — > E — hm only for the activated rates However, it has been advocated also as an approximate 
treatment for e.g. Glauber rates © J0|, and so is worth considering for general Wq(E' *— E). 
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From ([13)1 . the first term in square brackets is seen to give —/3(dC(t,t w )/dt. For the 
second one, one notes from (j!3|) and (0) that 

dC 



dt v 



•r (-B)(t-t w ) 



T (E)P (E,t w ) + -^-(E,t y 



(24) 



JdEA 2 (E)e~ 

= jdEA 2 (E)e- Vomt - t ^p(E) JdE'w (E <- E')P (E',t w ) (25) 

which apart from pref actors is just the second term in (f2*3|) . Thus, for an?/ mean-field trap 
model with the multiplicatively perturbed rates (|T§|) . and any zero- mean observable, one 
has the result given by Bouchaud and Dean jlj for Bouchaud's trap model 

R(t,t w ) = -P(^+P(l-0^- (26) 
ot at w 

The above calculation shows that this relation holds entirely independently of the precise 
form of the trap depth distribution§ p(E) or the transition rates Wq(E' <— E). In 
equilibrium, where C(t,t w ) is a function of t — t w only, it of course recovers the usual 
FDT, R(t,t w ) = (3 dC(t,t w )/dt w . Equation (ffilj) applies in particular to (zero-mean) 
neutral observables [HJ dj, where m is completely decoupled from E and therefore 
p(m\E) is independent of E. It remains true also for more general observables, however, 
as long as they have zero conditional mean ©. 

5. The Barrat and Mezard model 

Next let us turn to the Barrat and Mezard model, with the natural prescription which 
assumes that the field shifts all energies. From (jSJ) the rates are then 

w(E' m' <- E,m) = -— — — — — - — -pr (27) 

v ' ' ' 1 + exp{f3[(E' - hm') - (E - hm)}} v ' 

For low T, Ritort 10 argued that as a reasonable approximation to this one could 
consider multiplicatively perturbed rates 

w(E', m! <-E,m) = e h ^ m '-^ m) w G {E' <- E) (28) 

with 7 = ji = 1/Tg for exponential p(E). Equation (|2^j) is identical to (fT^j) apart from 
the replacements — () — > 7, (3C, — > p,. As expected from the general result (|27)j) for 
multiplicatively perturbed rates, Ritort therefore obtained the relation 

R {t , tw) = - li - +7 - (29) 

between response and correlation. This was found confirmed in simulations. However, 
as discussed in the appendix, these simulations were effectively performed directly with 
the approximate rates ()28|). so did not give a check of how well this approximation 
captures the behaviour of the Barrat and Mezard model. I now show that the response 

§ The irrelevance of the form of p(E) may well have been known to the authors of Ref. 0], but was 
not stated there. The version of (|26|l given in 0] is nevertheless more limited than the one given here, 
since only activated rates and neutral observables m were considered. 
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Figure 1. Effect of a field on Glauber transition rates, sketched for h(m! — m) > 0. 
Solid lines show the original transition rates (JSJ, dashed lines those in the presence 
of a field h, which are shifted to the right by h(m' — to); see arrow on the left. The 
difference between the two curves is the change in the rates, equation i|30|l : its integral 
is clearly h(m' — to). Left: Case where \h(m' — m)\ > T; the range where the change 
is significant is given by h(m' — to). Right: Case where \h(m' — m)\ < T; here the 
temperature T sets the range where rates change significantly. For small h and small 
T the range is small in either case. 



can be calculated exactly even with the exact rates (|27|). and that the results differ from 
those found for multiplicatively perturbed rates. 

To calculate the response function, consider the change in the transition rates due 
to the field, 

Aw(E', w! <-E,m) = - - exp{(3[{E , _ w) _ {E _ hm)]} ~ i + exp[/3(jE; , _ E)] ( 30 ) 

This is significantly different from zero only for E' within a range of order max{T, \ h(m'— 
m)\} around E; see figure [D If this range is small compared to T g , which is true for 
T T g and small fields h, then in 

AT(E, m) = j dE' dm p(m'\E')p(E')Aw(E' , w! «- E, m) (31) 

we can to leading order replace E' by E in the factor p(E'); the same is true for the first 
factor if we assume that p{m'\E') varies with E' at most on the same scale (~ T g ) as 
p(E'). Using / dE' Aw(E', m! E, m) = h(m' — m), which from figure^is geometrically 
obvious, together with (0) one thus finds 

AT(E, m) = J dm' h(m - m)p(m'\E)p(E) = -hmp(E) (32) 

The same argument can be applied to the integral in the second term of ([18]). as long 
as we are in an out-of-equilibrium regime where Po(E',t w ) varies with E' on a scale of 
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T g , rather than T as it would in equilibrium. This gives to leading order 
JdE' dm' Aw(E,m <- E', m')p(m'\E')P (E' , t w ) = 

= J dm' h(m - m')p{m'\E)P {E, t w ) = hmP (E, t w ) (33) 

One can now insert ()32I33|) into (|18|): after carrying out the m-integration and 
simplifying one sees that both terms give the same contribution. Dividing by h, the 
linear response function is therefore 

R{t, t w ) = 2 J dE A 2 (E)e- r ° iE)it - t - ) p(E)P (E, t w ) (34) 

Although I had implicitly assumed an exponential p(E) above, this result obviously 
remains valid also for other p(E), as long as T is much smaller than the energy scale 
over which p(E) and p(m\E) vary significantly. Comparing with (|13|). one now sees that 
in general there is no simple relation between the response and correlation functions 
for the Barrat and Mezard model. However, for the exponential trap distribution 
p(E) = T g _1 exp(E/T g ) such a relation does exist. For low T one can approximate 
the transition rates by a step function, wq(E' <— E) « 0(i? — E') and the total exit 
rates are 

T (E) = f E dE'p(E') = e E ^ = T g p(E) (35) 

J — oo 

Thus, comparing ffTTTft and (J33|) gives 

9 8(1 

R ^—T g m (36) 

Surprisingly, this is not dissimilar to the result (|26|) which one obtains for Bouchaud's 
model in the most natural case £ = 1: the only difference is in the prefactor, which is 
1/T for Bouchaud's model but 2/T g for the Barrat and Mezard model considered here. 

A simple application of (j3*Uj) is to the case of a neutral observable, with A 2 (E) = 1 
(say) independently of E. Then from (fT3|) one sees that C(t,t w ) is the hopping 
correlation function, i.e. the probability of not leaving the current trap between t w 
and t. This was worked out by Barrat and Mezard [2] for T — ► 0, with the result that 
C(t,t w ) = t w /t for long times. Equation (jSfiJl then yields R(t,t w ) = (2/T g )t w /t 2 ; the 
step response follows as 



x(Mw) = f ' dt' R(t,t' 



* 1 

-L p- 



1 



T 



l-C 2 (Mw) (37) 



An FD plot of x versus C therefore has a parabolic shape, with vanishing asymptotic 
slope dx/dC for C — > 0, i.e. = 0. The above calculation shows that the result (|3*7j) 
is exact for the Barrat and Mezard model in the limit T —>■ 0; it is also consistent with 
simulation results as shown in figure El 

We can now assess the accuracy of the approximation of multiplicatively perturbed 
rates (|2*Hj) . From (f2"H|) . one finds in this case (THj, by arguments analogous to those 
above, that x — 7(1 — C) + f (1 — C 2 ). Recalling that 7 = p = 1/T g , this is seen 
to be rather different from (J37|) . In particular, the approximation of multiplicatively 




Figure 2. Correlation and response of neutral observables, for the Barrat and Mezard 
model at T = 0; energies are scaled such that T g = 1. Main plot: C(t, i w ) and x(t, t w ) 
against t, for t w = 50. The response was determined for field h = 0.1, which can be 
checked to be in the linear regime. Circles: simulation results (sec appendix). Lines: 
theoretical predictions C — t w /t and x = 1 — {tw/t) 2 , see equation (|S7)l . Inset: FD 
plot of x vs C- 



perturbed rates incorrectly predicts a non- vanishing asymptotic slope of the FDT plot, 
dx/dC = - 7 = -1/T g . 

Finally, it is worth discussing a difference between the T — > correlation and 
response functions for neutral observables, in terms of their dependence of p(E). The 
hopping correlation function is independent of p(E), as shown in [2]. Within the 
formalism used here, this is clear if in (jlHj) one sets A 2 (E) = 1 and changes variables to 
the cumulative trap density r(E) = ff^ dE' p(E'). Together with the first part of 
this gives C(t,t w ) = J^dr e _r ^ _ * w ' ) Po( r , t w ). Since Po(r, t w ) is independent of p(E) (as 
can be shown from (J2J) using the same change of variable), the same then holds for 
C(t,t w ). The intuitive reason for this independence is that the T — > Glauber rates 
wq{E' <— E) = Q(E — £") depend only on the relative "height" of departure and arrival 
trap, but not otherwise on the actual values of E and E'\ correspondingly, the total exit 
rate Tq(E) depends only on how many traps are at energies below E, i.e. on r(E). 

By contrast, the response function does depend on p(E): transforming from E to 
r in gives R(t,t w ) = 2 J^dr e~ r ^ iw ) p(E(r))P (r,t w ) and the dependence on p(E) 
cannot be eliminated. This can be explained intuitively by noting that the perturbation 
term —hm which shifts the energies E introduces an energy scale which is not present 
for h = 0. The response is sensitive to how many traps there are with energies near 
(measured on this scale) that of the departure trap, and hence to p(E). 
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6. Conclusion 

In this paper I have considered mean-field trap models, which are simple and intuitive 
models of glassy dynamics. I showed briefly that a relation between out-of-equilibrium 
correlation and response functions in these models, first given by Bouchaud and Dean, 
is valid for a general class of mean-field trap models; it requires only that the transition 
rates are affected in the simple multiplicative way (fT5j) by an applied field. 

I then considered the Barrat and Mezard model, which has Glauber dynamics and 
an exponential distribution of trap depths. Glassiness arises in this model from entropic 
barriers, rather than energetic ones as in the Bouchaud model, and so it is of interest 
to compare the FD relations that result from these different physical mechanisms. In 
the natural version of the model where the effect of a field is to shift the energies of 
all traps according to the usual prescription E — > E — hm, the effect on the transition 
rates is not simply multiplicative. The out-of-equilibrium response can nevertheless 
be obtained exactly for low T, and one finds a relation which is quite similar to, but 
distinct from, that given by Bouchaud and Dean. The exact calculation also shows that 
an approximate treatment using multiplicatively perturbed rates [10] gives qualitatively 
incorrect results. 

Comparing the above results for the (natural) Barrat and Mezard model with those 
for Bouchaud's model (with, likewise, the natural choice ( = 1), one notes two intriguing 
parallels for the low-temperature out-of-equilibrium dynamics. Firstly, both models 
give FD plots with = 0, i.e. with a slope dx/dC which tends to zero in the limit 
C — » 0. (For the Barrat and Mezard model with non-neutral observables|| this can be 
deduced by applying the arguments of [TU] to the relation Second, the value of the 

susceptibility itself in the same limit is Xoo = i n both models for neutral observables; 
see again (|57j). This is precisely the value that one would expect if, as T is lowered, 
Xoo "freezes" at T = T g and remains independent of T for T < T g . For Bouchaud's 
model this T-independence can indeed be shown JU]; for the Barrat and Mezard model 
the result Xoo = f° un d above for T — > strongly suggests that Xoo is likewise 

T-independent for < T < T g . Even though the slow out-of-equilibrium dynamics in 
the two models is very different, being caused by activation over energy barriers for 
Bouchaud's model and by entropic barriers for the Barrat and Mezard model, we thus 
have the intriguing observation that some features of the out-of-equilibrium FD relations 
are shared. It will be interesting to explore whether this correspondence extends to other 
properties, and possibly to other models of glassy dynamics. 

|| Strictly speaking, as shown in |E3], one requires observables that probe only the aging dynamics, in 
the sense that their correlation function only decays on timescales that grow with t w . A counterexample 
would be an observable that is sensitive only to the very shallow traps, which in Bouchaud's model 
gives a correlation function that decays completely on timescales of 0(1) 0. 
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Appendix: Simulation method 

To simulate any mean-field trap model, one can use that in the limit iV — > oo no trap is 
visited twice, so that E and m can be sampled anew at each transition and no explicit 
population of traps needs to be maintained. The probability for making a transition 
from a trap with (E,m) to one with (E',m') is 

P 1 (E',rri\E,m) = T- 1 ( y E,m)p(m'\E')p(E')w( y E',m' <- E,m) (38) 

and contains T(E, m), the total exit rate from the current trap, as a normalization factor; 
see (JSJ). Now specialize to the Barrat and Mezard model, with p(E) = T" 1 exp(£ , /T g ), 
E < and a neutral observable for which I take p(m\E) = p(m) as a zero mean, 
unit variance Gaussian independently of E. The transition rates at T = are 
w(E', ml <— E,m) = Q(E — hm — E' + hm'). Integrating over E' in (jHSj) then gives 

P x {m'\E y m) = T'\E, m)p{m')e {E - hm+hm ' )/T ^ oc p{m')e hm ' (39) 

Dividing by this yields 

P x {E'\ml, E, m) = G(E - hm + hm' - E ')T- 1 e {E '- E+hm - hm ' )/T ^ (40) 

which is just an exponential distribution over — oo < E' < E — hm + hm'. One can thus 
sample from ()38|) by first sampling ml from (|39|) , which is a Gaussian with unit variance 
and mean h; after that one samples E' from (|4()jl . The total exit rate follows e.g. from 
normalization of (}3*9*j) as 

T(E, m) = exp[(E - hm)/T g + h 2 /(2T*)} (41) 

It is important to note from (|4()jl that the distribution of E' depends on m — 
m! . One might be tempted to neglect this dependence, replacing (}4T)j) by Q(E — 
E')T- 1 e( E '- E)/T * [Ej. However, by repeating the calculations leading to (j39l4()l41jl 
one easily checks that this is equivalent to changing from the exact rates (|27j) to the 
multiplicatively perturbed rates with 7 = p = 1/T g . As shown above, this leads to 
rather different response functions; the precise form of (}4*0|) is thus important to get the 
correct results. 

The results shown in figure El were obtained for a quench from T = 00 at t = 0, 
corresponding to the initial condition P(E, m, 0) = p(m)p(E), and averaged over 5 x 10 7 
runs. Direct simulations with a population of iV = 10 8 traps yielded equivalent results, 
though one needs to be aware of finite- iV effects which become more acute for low E 
because of the exponential decrease in the population density p(E). 
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